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Abstract 

We construct an effective field theory valid for processes in which highly 
energetic light-like particles interact with collinear and soft degrees of freedom, 
using the decay B — » X s j near the endpoint of the photon spectrum, x = 
2Ery/mb — ► 1, as an example. Below the scale [i = mj both soft and collinear 
degrees of freedom are included in the effective theory, while below the scale 
[i = rriby/x — y, where 1 — y is the lightcone momentum fraction of the b 
quark in the B meson, we match onto a theory of bilocal operators. We show 
that at one loop large logarithms cancel in the matching conditions, and that 
we recover the well known renormalization group equations that sum leading 
Sudakov logarithms. 
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I. INTRODUCTION 



Effective field theories (EFT's) provide a simple and elegant method for calculating 
processes with several relevant energy scales Part of the utility of EFT's is that they 
dramatically simplify the summation of powers of logarithms of ratios of mass scales, which 
would otherwise make perturbation theory poorly behaved. For example, in a theory with a 
very heavy particle of mass M, one-loop corrections will typically be enhanced by log(M/ A), 
where A is a low scale in the problem. In the EFT in which the heavy particle has been 
removed from the theory, such logarithms are replaced by factors of log(/z/ A) (where n is the 
renormalization scale in dimensional regularization, or the cutoff in cutoff regularization), 
and the complete series of leading logarithms a™ log n (/i/ A) is straightforward to sum via the 
renormalization group. 

The situation is more complicated for processes with highly energetic light particles. In 
this case, there are both collinear and infrared divergences in the theory, which give rise to 
the familiar Sudakov double logarithms 0. For example, the perturbative expansion of the 
iV'th moment of the photon spectrum in inclusive b — > X s 7 decay is of the form 

E E C n , m a:\og m N. (1) 

n m<2n 

Although the arguments of these logarithms are not obviously the ratio of two scales, they 
arise because the typical energy and invariant mass of light particles are widely separated, 
and they may be summed via well-known techniques based on factorization theorems |3j 
into the form 



exp 



E (a n a n s log n+1 N + b n a n s \og n N) + 



(2) 



The terms <log n+1 iV are referred to as the leading logarithmic contribution, the terms 
a™ log n iV are referred to as the next-to-leading logarithmic contribution, and the remaining 
terms are called subdominant. 

Recently there has been some discussion in the literature of summing Sudakov loga- 
rithms using effective field theory techniques Such an approach could have several 
advantages over the conventional method; in particular, while factorization formulas are 
based on perturbation theory, EFTs, by construction, are valid beyond perturbation theory, 
and by including higher dimension operators it should be straightforward (if tedious) to go 
beyond the leading twist approximation. In the various versions of the EFT approach which 
have been suggested, the effective theory is the so-called "Large Energy Effective Theory" 
(LEET) which describes light-like particles coupled to soft degrees of freedom. However, 
a difficulty with the approaches presented to date is that, as pointed out in Refs. 0, in the 
minimal subtraction (MS) scheme logarithms arising at one loop in LEET do not match 
logarithms arising at one loop in QCD for any choice of the matching scale fi; hence these 
logarithms may not be summed using the RGE's. 

In this paper we consider this problem in the context of B — > X s 7 decays.0 We show that 
the problem of matching scales may be resolved by introducing a new intermediate effective 



In fact, the authors of || argued that the resummation of subleading Sudakov logarithms is not 
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theory containing both soft and collinear degrees of freedom, which is then matched onto 
LEET (effectively integrating out the collinear modes) at a lower scale. We show that the 
matching conditions onto both effective theories contain no large logarithms at one loop. 
We then calculate the RGE's in the two theories summing the leading logarithms and a 
certain subset of the next-to-leading logarithms. To this order the expression obtained for 



the resummed Sudakov logarithms is identical to that derived in Refs HlO 



II. SUDAKOV LOGARITHMS IN B -> X s -f AND LEET 

Inclusive decays of heavy quarks have been well understood for many years in the context 



of an operator product expansion (OPE) in the inverse mass of the heavy quark [11]. At 
leading order in the Aqcd/^6 expansion the B meson decay rate is equal to the b quark 
decay rate, and nonperturbative effects are suppressed by at least two powers of AqcdA?V 
However, the OPE only converges for sufficiently inclusive observables. Unfortunately, ex- 
perimental cuts on measurements of rare decays such as B — > X s 7, B — > X s £ + £~, and 
B — > X u £u are required, restricting the available phase space considerably. Since all of these 
decays are of phenomenological interest, either in the determination of \V u b\ or detection of 
new physics, understanding inclusive decays in restricted regions of phase space is important. 

If the phase space is restricted such that the final hadronic state is dominated by only a 
few resonances, the breakdown of the OPE simply reflects the fact that an inclusive treatment 
based on local duality is no longer appropriate. This is the case for the dilepton invariant 



mass spectrum in inclusive B — ► X s £ + £~ and B — > X u £v decays ||12|| . However, when the 
kinematic cut is in a region of phase space dominated by highly energetic, low invariant 
mass final states, the OPE breaks down even for quantities smeared over a parametrically 
larger region of phase space, where the decay is not resonance dominated. This situation 
arises in the endpoint region of the electron energy spectrum and the low hadronic invariant 
mass region in semileptonic B — > X u £v decay, as well as the endpoint region of the photon 
spectrum in B — > X s ^ decay |PJT3[] . 

Consider the dominant contribution to the decay B — > X s 7, which arises from the 



magnetic penguin operator |14 



d > = 7T^ m <> s^t 1 + 75)6 F, u , (3) 

l07T Z 1 

where the strange quark mass has been set to zero.0 The OPE for this decay is illustrated 
in Fig. We write the momenta of the b quark, photon, and light s quark jet as 

f h = m bV » + <f = yin", p» = + P + F (4) 

where, in the rest frame of the B meson, 



necessary for practical purposes for this decay. Nevertheless, it provides a simple example in which 
we may compare our results to those in the literature. 

throughout this work we will ignore the contribution of operators other than O7 to the decay. 
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1,0), n" = (1,0, 0,-1) 



rr 



[1, 0,0,1). 



(5) 



Here k^ 1 is a residual momentum of order Aqcd, and P = -^(1 — where x = 2E 1 /m b . 

The invariant mass of the light s-quark jet 



p 2 ~ m b n ■ (I + k) = ml(l — x + 



(6) 



(where &;+ = k + /m b ) is O(m^) except near the endpoint of the photon energy spectrum 
where x — > 1. Inclusive quantities are calculated via the OPE by taking the imaginary part 
of the graphs in Fig. [I] and expanding in powers of k^/ ytf%- As long as x is not too close to 
the endpoint, this is an expansion in powers in k^/mt, which matches onto local operators. 
This leads to an expansion for the photon energy spectrum as a function of x in powers of 




+ . 



P b =m b v+k 



a s and l/m b |L5] 



FIG. 1. The OPE for B -» X s7 . 
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and the subscript "+" denotes the usual plus distribution, 



1 



(1-x) 
' log(l -x) \ 

, (1-*) ) 



lim , 

13^0 I 1 — X 

lim( lQg(1 ~ X V l - * - /3) + ~ hg 2 (P)5(l -x-0) 

13^0 1 — X 2 



(7) 



(8) 



(9) 



The parameters Ai and X 2 are matrix elements of local dimension five operators. 

Near the endpoint of the photon spectrum, x — > 1, both the perturbative and nonpertur- 
bative corrections are singular and the OPE breaks down. The severity of the breakdown 
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is most easily seen by integrating the spectrum over a region 1 — A < x < 1. When 
A < ^QCT>/ m b the most singular terms in the 1/rrif, expansion sum up into a nonperturba- 
tive shape function of characteristic width AQcn/ m b ||16|| . The perturbative series is of the 
form 



1 f 1 dT 

r A-a dx 



ot K C 



1 + ^ (-21og 2 A - 71ogA + ...)+ O (af) 



(10) 



where the ellipses denote terms that are finite as A — ► 0. These Sudakov logarithms are 
large for A < 1, and can spoil the convergence of perturbation theory. The full series has 
been shown to exponentiate [||[10| and the leading and next-to-leading logarithms must be 
resummed for A < exp (— yV / a s (rrib)^ , which is parametrically larger than A QCD / m h in the 
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In general, "phase space" logarithms are to be expected whenever a decay depends on 
several distinct scales. For example, in b — > X c ev e decay the rate calculated with the OPE 
performed at fi = m b contains large logarithms of m c /mfc. In |1| an EFT was used to run 
from rrif, to m c , summing phase space logarithms of the ratio m c /m(,. Similarly, in b — > X s 7 
near the endpoint of the photon energy spectrum the invariant mass of the light quark 
jet scales as m^l — x, and is widely separated from the scale /i = m b where the OPE is 
performed. In order to sum logarithms of A (or the more complicated plus distributions 
in the differential spectrum, Eq. (|7|)) we would expect to have to switch to a new effective 
theory at fi = rrib, use the renormalization group to run down to a scale of order m^X — x, 
at which point the OPE is performed. (In fact, we will see that the situation is slightly more 
complicated than this). 

We are then left with the question of the appropriate theory below the scale m^. The 
simplest possibility is to expand the theory in powers of k^/rrib and l^/rrib. The heavy 
quark is then treated in the heavy quark effective theory (HQET) ||19H , while the light quark 
propagator is treated in the large energy effective theory (LEET) proposed many years ago 
by Dugan and Grinstein |7]]. Expanding the s quark propagator in powers of l/m b) we find 
the LEET propagator 



p 2 s 2 n ■ (I + k) 



O 



. m b , 



(11) 



LEET is an effective theory of lightlike Wilson lines, much as HQET is an effective theory 
of timelike Wilson lines 0]. The hope would then be to match QCD onto LEET and then 
use the renormalization group to sum the Sudakov logarithms. This is the approach taken 
in . However, a simple attempt at matching shows that this does not sum the appropriate 
logarithms. 

Consider the one-loop matching of the operator Oj from QCD to LEET. We regulate 
ultraviolet (UV) divergences with dimensional regularization (d = 4 — 2e). We introduce a 
small invariant mass p 2 s for the s quark which regulates all infrared (IR) divergences except 
that in the heavy-quark wave function diagram, Fig. 0(b). This IR divergence is regulated 
using dimensional regularization. The vertex diagram, Fig. 0(a) yields 



= -CMsT»b 



47T 



log 2 4 + 21og4 + --- 
m b m b 



(12) 
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where 



8tt 



: m b cr 



(1+75) 



-<lu 



(13) 



C^dj,) is the Wilson coefficient of O7 and the dots denote (here and in the rest of the paper) 
finite terms which are not logarithmically enhanced. Including a factor of \fZ for each 
external field 



1 - 
1 



a,C. 



s^F 

Att 
ol s Cf 

47T 



- + 31og4 + 

e mi 



1 v 2 u, 

--log^_+logA r + ... 



mi 



mi 



where 

ft 2 = 47r/iV 7B 

and adding the counterterm required to subtract off the UV divergence 

Ot s CF 1 



1 + 



An 6 



we find 

Aqcd = C 7 (fi)sT% 



a s C F (. 2 pi 3 , p 2 1 _ . u 2 
' log 2 ^ + - log ^ + - + 2 log " 



An V ° m 2 ' 2 



mi 



nn. 



(14) 
(15) 

(16) 

(17) 

(18) 




FIG. 2. One-loop corrections to the matrix element of O7 in QCD. 

The corresponding LEET diagram is shown in Fig. || Neither of the wave function graphs 
gives a contribution, since the light quark wave function in Feynman gaugeQ is proportional 
to n 2 = 0, and the heavy quark wave function vanishes in dimensional regularization. Thus 
the only contribution is from the vertex graph. Denoting the coefficient of the corresponding 
operator in LEET as C^(l + (a s C F /An)C^ + ...), we find 



3 We will work in Feynman gauge throughout this paper. 
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FIG. 3. One-loop correction to the bs^f vertex in LEET. 



where p 2 /rrib is the soft scale. So 




Li 2 m 2 b \ T(l + 2e)r(l - 2e)r(l + e) 



ft 



2, ft 
-log r 



2 „2 

21og 2 ^ + . 



(19) 



C7(°)(/2) = C7 7 ( / i) 



rribfx 



0+21og 2 ^l-log 2 4 



3, ft 2 
-log-^ 



-2 log 4 + 



(20) 



We immediately notice two problem^: 



1. There is no matching scale fx at which all the large single and double logarithms in 

vanish. Thus, there are logarithms in the rate which cannot be summed using the 
renormalization group in LEET. 

2. C« contains a divergence proportional to Mogp 2 . Since p 2 is an infrared scale in 
the problem, it is not clear how to sensibly renormalize this term. In Ref. this 
divergence was cancelled by a nonlocal counterterm in the inclusive rate; however, this 
term indicates that LEET cannot be used for exclusive processes p0 |. Furthermore, 
the matching of the inclusive rate performed in || still leaves large logarithms in the 
coefficient of the operator. 

The problem is that LEET only describes the coupling of light-like particles to soft 
gluons, but does not describe the splitting of an energetic particle into two almost collinear 
particles. Thus, by matching onto LEET one is integrating out the collinear modes which 
also contribute to infrared physics. As we will show below, once collinear degrees of freedom 
are included, both of the above problems are resolved. 



4 Note that CV(/i) includes a factor of a s Ci?/(47r) log(mw/Li), which converts one of the factors of 
log(/«/m&) in Eq. (|l^) to log(mw / 'mb) ■ This is not important for our argument. 
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III. THE COLLINEAR-SOFT THEORY 



A. Collinear and soft modes 

It is convenient to work in light-cone coordinates p^ = (p + , p~ , p\) , where p + = n-p and 
p~ = n-p, and to define a power-counting parameter A = y/1 — x that becomes small in the 
limit x —>■ 1. The momentum of the light-quark jet then scales as 

^~m b (A 2 ,l,A). (21) 

This scaling is unchanged by emission of either soft or collinear degrees of freedom, with 
momenta scaling as 

Psoft ~ m b (X 2 , A 2 , A 2 ), Pcoiiinear ~ m b (X 2 , 1, A), (22) 

and so emission of both modes is kinematically allowed. It is the presence of infrared sensitive 
graphs with collinear loop momentum that makes this EFT more complicated than other, 
more familiar, EFT's, where infrared sensitivity comes purely from soft modes. This is 
similar to the situation in non-relativistic QCD (NRQCD) in which power counting is 
complicated by the fact that a given amplitude receives contributions from loop momenta 
which are small compared to the heavy quark mass, but which have parametrically different 
dependence on the heavy quark velocity v. In NRQCD, the relevant scales are known as 
soft, ultrasoft and potential, and must be treated separately in order to obtain consistent 
power counting p2| , |23| . 

We follow a similar approach here, and introduce separate fields for both soft and collinear 
degrees of freedom^. Between the scales m& and m b X the effective theory contains separate 
fields for both collinear and soft modes, while at scales below ~ m b X (the exact scale depends 
on the operator under consideration, as will be discussed in the next section), the collinear 
modes are integrated out of the theory and it is matched onto LEET. We will refer to 
this intermediate theory as the collinear-soft theory, and resist the urge to create another 
acronym. 

There is an important difference between the approach taken here and the one taken in 
Refs. P^ , |25| where logarithms of v are summed in NRQCD and NRQED. In the latter case 
no intermediate theory is introduced; instead the running is performed in one step through 
the velocity RGE. In NRQED these two approaches differ at subleading order [25,26|, and 
it may be that such one-step running is needed here at two loops. 

The power counting rules in the collinear-soft theory may be obtained by a field rescaling, 
analogous to that performed in p^| . The scaling of the fields is chosen such that the 
propagators are all 0(1), putting the A dependence into the interaction terms. For example, 
in the kinetic term for a soft gluon, 

d 4 x(d,A a u -d„A a \ 2 (23) 



5 At two loops an additional gluon field scaling as (A, A, A) might have to be included 
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Factor 


Scaling 


soft gluon 


A 2 


collinear gluon 


A 


heavy quark h 


A 3 


collinear quark £ 


A 


collinear volume element d^x c 


A- 4 


soft volume element d x s 


A- 8 



TABLE I. Power counting rules for fields in the collinear-soft theory in Feynman gauge, where 



X = VT 



l 2~pT > collinear quark ~ijT ^IRRRP collinear gluon 

iP 2^ v 

> heavy quark ~i TT HTfRRP soft gluon 



v-p 



FIG. 4. Propagators in the collinear-soft effective theory. 



the typical length scale associated with soft excitations scales as A~ 2 ~ p^,ft> so the factor 
of d A x scales as A -8 . Each derivative scales as p so ft ~ A 2 , so the soft gluon field must scale 
as A 2 for the kinetic term to be 0(1). 

Since the various collinear momentum components scale differently with A, power count- 
ing for collinear gluons is gauge dependent (this is easily seen from the propagator, since 
in a covariant gauge the components of the k^k v term scale differently). In this paper we 
are working in Feynman gauge, in which case the different components of collinear gluons 
have the same scaling. Performing a similar analysis for the other fields, we obtain the 
power-counting rules given in Table |I[ 

Rather than write down the effective Lagrangian for the various fields, which is quite 
lengthy, we will instead just give the Feynman rules, which are obtained by expanding the 
QCD amplitudes in powers of A. The spinors in the collinear-soft theory are related, at 
leading order in A, to the QCD spinors via 

h v = P+u , £ n = P n u , i n = PfiU , (24) 

where we have defined the projection operators 

p_^ +1 p _ # p _ # (n^ 

r + — 2 ' n — X " — ~4~ ' 

which project out the heavy quark spinor, a massless spinor in the n direction, and a massless 
spinor in the n direction respectively. The propagators for the different fields are shown in 
Fig. § 

The interactions leading in A which we will need in this paper are shown along with 
their Feynman rules and scaling in Fig. [5| Note that the interaction of a soft particle with 
a collinear particle leaves the minus and perpendicular momenta of the collinear particle 
unchanged, since they are parametrically larger for the collinear particle. This is analogous 
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(a) (b) (c) 

FIG. 5. Leading order quark-gluon interactions in the collinear-soft effective theory: (a) 
collinear-collinear, (b) collinear-soft, and (c) heavy-soft. Applying the rules from Table |, the 
vertices scale as (a) A -1 , (b) A , and (c) A . 



to the multipole expansion which is performed in NRQCD ||281 . As a result at one loop 



soft-collinear interactions in this theory are equivalent to LEET, since collinear propagators 
in soft loops reduce to LEET propagators: 

ft n-(p-k) ft pr = _\__ 

2 (p-k) 2 ~ 2 (p - k)+p- - (p 1 ) 2 2n-k 1 ' 

where p is a collinear momentum, k is a soft momentum, and p 2 = from the equations of 
motion. Once again, this is analogous to NRQCD, where in ultrasoft loops the Feynman rules 
reduce to those for HQET. By the same token, in soft-collinear interactions, the appropriate 
volume element is the collinear volume element, scaling as A~ 4 . 

Because the leading purely collinear interaction, Fig. |5|(a), scales as A -1 , power counting 
for collinear loops is less simple than for soft loops. Terms which would scale as A -2 , such as 
the purely collinear wave-function graph in Fig. || are proportional to n 2 = and so vanish 
in the effective theory. However, the 1/A coupling enhances terms which would naively be 




FIG. 6. Pure collinear wavefunction graph. The 0(X 2 ) contribution vanishes. 

suppressed. In fact, although the A-counting looks complicated, graphs with only collinear 
lines are identical to the corresponding graphs in QCD. This is because in any graph in 
which all the lines have the same scaling (and there are no purely soft graphs, so this only 
refers to purely collinear graphs), expanding in powers of A does not change the propagators. 
Since the locations of poles in the propagators are unaffected, it is irrelevant whether one 
calculates the full graph in QCD and then expands in powers of A, or calculates each order 
in A in the collinear-soft theory. Thus, for purely collinear graphs, such as the wave function 
graph in Fig. |6|, we will not bother to write down the complete set of operators, but simply 
calculate the graph in QCD and expand. 

There is one important subleading operator, shown in Fig. ^, which can be enhanced 
by the 1/A piece of the purely collinear coupling. By momentum conservation, there is no 
vertex coupling two heavy quarks and a collinear quark, since a heavy quark cannot emit 
a collinear gluon and stay on its mass shell. However, expanding the diagram in Fig. [7] in 
powers of A gives the nonlocal 0(A) interaction shown in the figure. (This is similar to the 
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nonlocal operators found in p9|). Though it is formally subleading, in graphs such as Fig. 
H(a) it gives an 0(1) effect. 




FIG. 7. Non-local vertex at 0(A). 



B. Matching onto the Collinear-Soft Theory 

We now proceed to compute the matching conditions for the operator 7 , and demon- 
strate that there are no large logarithms in the matching coefficients. At tree level, the 
matching is trivial. Defining the current in the effective theory by 

W = l n T»h v , (27) 



where T M is given in (|13"D , the Wilson coefficient Cy at tree level is 

C v = l + 0(a s ). (28) 

To perform this matching at one-loop, we repeat the one-loop matching calculation discussed 
in Section [HJ but now using the collinear-soft theory instead of LEET, hence including 
collinear modes. The calculation is simplest if we set the invariant mass of the s quark to 
zero; this introduces additional infrared divergences to the calculation which cancel in the 
matching conditions. The one loop matrix element of O7 in full QCD can be calculated from 
the diagrams in Fig. ^ and we find the amplitude 



^4qcd = sr"6 



_ n log(A 2 /mg) 5_ 1 , 7 l_ 

4tt L 2 + e + 2e 2 g m 2 , 2 g ml 



(29) 



where all the 1/e divergences are infrared in origin. The one loop correction in the collinear- 
soft theory can be calculated from the Feynman diagrams in Figs. ^| and || In pure dimen- 
sional regularization all graphs are zero, as there is no scale present in the loop integrals. 
Thus, we find the matching condition 



CyZy = l + ^ ±+ ^ ^ + f + ■■■ , (30) 




where Z v is the counterterm required to subtract the UV divergences in the collinear-soft 
theory. 

This derivation of course assumes that the collinear-soft theory reproduces the infrared 
behaviour of QCD. We can check this by instead introducing a small invariant mass for the 
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s quark, as in Section ||, and explicitly verifying that the dependence on the invariant mass 
in the collinear-soft theory is identical to that in full QCD given in Eq. fll8l) . The soft gluon 
contribution in the collinear-soft theory is identical to the LEET result, given in ([19]) 




= Cvl n T»h- 



-C v £ n T»h 





2n \ 


pi 


a s C F 


2 


An 


~7 



r(i + e)r(l - e)r(2 - e) 1 



T(2 - 2e) 



2 

e 



A 2 



log 2 ?| + 2 log 4 + 



(32) 



As previously discussed, the leading piece of the wave function graph Fig. |](b) is 0(1/A 2 ), 
but fortunately vanishes. Expanding to higher orders in A, the graph gives the same result 
as in full QCD, ([15]). We therefore obtain for the contribution of the collinear gluons 



An 



2 3 2, p\ 
2e e n z 



i 2 Pi 3 1 Pi 
log 2 ^ + -log 



Adding the soft and collinear contributions, as well as the counterterm given in 
obtain 



A r 



-C v £ n T»h 



a s C[ 
An 



i 2 Pi 3 . pi 1 1 , 2 V 2 3 , ^ 2 

log 2 i^ + -log^ + -- - log 2 ^ - - bg + . . . 

mi 2 mi e 2 mi 2 mi 



(33) 



we 



(34) 



Note that the troublesome divergence ~ -\ogp 2 s cancels once the two contributions 
and (|33|) are added. Thus, both collinear and soft modes are required for the theory to be 
renormalized sensibly. Comparing to the full theory result (0), we see that the collinear- 
soft theory reproduces the IR physics of QCD, and that at the scale p, = m b all nonanalytic 
terms vanish. This determines the matching scale to be mj, confirming the result found by 
calculating in pure dimensional regularization (BO). 



C. Renormalization group equations 

From the counterterm given in (|30|) it is simple to extract the anomalous dimension of 
the operator W* in the collinear-soft theory. From the definition, 
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lr = Z ?Ll_ + ^\z v (35) 



we have 



~®_Z = a s(^)C F 
dp, V 2ne 

d a s {P)C F (1 jl 2 5\ 



where we have used /3 = — ge + 0(g). This give the anomalous dimension 

a s (jl)C F ( jl 2 5\ 
^ = 2^( l0g ^! + 2j- (37) 

Note that the divergent piece of the anomalous dimension cancels between the two terms 
@. The RGE for the coefficient of the operator V 1 is therefore 

fij-Cvifr) = Tv{fi)CviP) ■ ( 38 ) 



Solving this RGE we obtain 
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C v tf) = [-^\ ^ C v (m b ), (39) 



where a s = a s (mh), flo = 11 — 2/3n/, and Cvijrib) = 1 + 0(a s (mb)). Note that in deriving 
the anomalous dimension (^) we have assumed that the nonlocal vertex given in Fig. |7| has 
the same running as the QCD coupling. This assumption needs to be checked in subsequent 
work. 



IV. THE SOFT THEORY 
A. Matching 



The collinear-soft effective theory is valid down to jl ~ m h \J\ — x, the typical invariant 
mass of the light s-quark jet. At this scale we integrate out the collinear modes, and perform 
an OPE to calculate the inclusive b decay rate. Diagrammatically, this is illustrated in Fig. 
(§). This results in a nonlocal OPE in which the two currents are separated along a light-like 
direction. As in Eq. (fj), we write the momentum of the eikonal line as 

P:=y^ + ^ + y(l-I/K (40) 

where is the residual momentum of the heavy quark (note that we distinguish y from x, 
the rescaled photon momentum, since beyond tree level they will differ). The imaginary piece 
of the first graph is then proportional to 5(1 — y + k + ) (where, as usual, hatted variables are 
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xj 




„ h(l-y-k + ) - 5(l-y-k + -q + ) 



FIG. 9. Diagrammatic representation of the OPE, as well as the zero and one gluon Feynman 
rules for the resulting nonlocal operator 0(y). 

divided by m^), so the OPE is in terms of an infinite number of nonlocal operators, labelled 
by y: 

0(y) = h v 8(l-y + iD + )h v . (41) 

Feynman rules for nonlocal operators of this type were obtained in |Kj , by writing them as 
the Fourier transform of operators in position space, and expanding out the path-ordered 
exponential in powers of the gauge field. Equivalently, the Feynman rules may be obtained 
by taking the imaginary piece of the time-ordered product in LEET with additional gluons; 
the single gluon Feynman rule is given in Fig. |9|. 

The matrix element of 0(y) between heavy quark states with residual momentum k is 

(b(k)\h v 5(l -y + ib + )h v \b{k)) = 5{l-y + k+) + 0(a s ) . (42) 



while its matrix element between hadrons is the well known structure function [16 



„ > {B\hJ(l - y + ,D + )h„\B) 

m = mMB) ' (43) 

Thus, LEET consists of a continuous set of operators labeled by y. Each operator has a 
coefficient that depends on the kinematic variable x, and the differential rate for B —>■ Xsj 
is given by the integral 

dT r 

— = r J dyC(y,x;/j l )f(y;n) , (44) 

where the C(y,x; /i)'s are the coefficients of the OPE. 

To match onto LEET at one loop we compare the differential decay rate in the parton 
model, b —>■ X s 7, which in LEET is 
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dr 

dx 



T dyC(y,x;fj)(b(k)\0(y;fj,)\b(k)). 



(45) 



We therefore need the one-loop matrix element of 0(y) between quark states. This may be 
calculated from the diagrams shown in Fig. [10. 




(a) 



(b) 



(c) 



FIG. 10. Feynman diagrams contributing to the one-loop matrix element of 0(y). 

Again all divergences are regulated in dimensional regularization. As an example, Fig. 
|T0l(a) gives 



(b{k)\0 {a \y)\b{k)) = iC F g 2 (JL\ 

\m b J 



4-d r 



dr 2 q± dq- dq + S(k + + 1 — y) — 5(k + + q + + 1 — y) 
(27r) d " 2 2tt 2tt (g + g_ - q\ + ie)(q+ + g_ + ie)q+ 

(46) 



The first term is proportional to 

d d q 1 
(2n) d (q 2 + ie) (q ■ v + ie)(q ■ n) 



1 rco 

dx I d\ 



A 



(47r) d / 



o Jo 



d d q 

~. / UX I CIA n 

(2n) d Jo Jo (q 2 + 2Xq-(v(l-x)+xn)y 

(47) 



ll r(3- d 2) j 1 dx f d\X< '((1 xf • ■>{[ ■. r A 



The A integral vanishes in dimensional regularization, so this term vanishes. After perform- 
ing the trivial q + integral in the second term, we are left with 



<&(A0|OM(l/)|6(*)>=i 



Cf9 2 



4-d 



1 



d d 2 q± dq. 



1 



2tt \m b J k + + 1 - y J (2n) d - 2 2n q_(k + + 1 - y) + q 2 ± - ie 

1 



C F g 2 / y, y- d 9Ck + + l-y) r 
2tt \m b J k,+l—y J 



q- - (k+ + 1 - y) + ie 
d d - 2 q ± 1 



y J {^) d - 2 ql + {k + + l-y) 2 
9(k + + l-y) 



mi 



(k + + i- y y+ 



2c 



(48) 



Using the identity 

6{y - x) 



(y-x) 



l+2e 



-—8(y-x)+B(y-x) 



(y-x). 



2e( l ^$) + 0(e 2 ) 



(49) 
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we find 

(b(k)\C)( a \y)\b(k)} — - 



a s C F 



47T { \ 

(1 — y + k 



1 1 , U 2 1 , 2 U 2 

-/ e mi 2 mi 



6(l-y + k + ) 



1 



m 



;i - y + k, 



/ log(l - 

, l-y + k+ 



(50) 



The diagram in Fig. 0(b) gives the same result as (a), while the diagram in Fig. |10|(c) 
gives 



(b(k)\0^(y)\b(k)) 



a s CF 

Alt 



?_ 21og 4), (1 _ s + i - + ) + 4^iM 
e m 2 ; + 



(51) 



In dimensional regularization the wavefunction diagrams vanish. Since the decay rate is 
infrared finite, including the wavefunction graphs simply converts an infrared 1/e divergence 
to an ultraviolet divergence. Therefore, we may neglect the wavefunction counterterm, and 
combining all graphs we find the bare matrix element 



(b(k)\O b ^(y) \b(k)) 



1 - 



a,C 



s^F 

4ir 



2 2 2, u 2 nl u 2 
--- + -log^-21og p 



mi 



mi 



, 2 
+ l0g 2 ^2 + 



5(1 - y + 



0(1 -y + fc H 



47T 



A 2 



4 log ^ + 4 



1 



(1-2/ + *+)+ 



+ 



'log(l - y + fc H 
, 1 - 2/ + K 



(52) 



where all divergences are ultraviolet. The renormalized operator 0(y;fi) is related to the 
bare operator by 



O h ^(y)= ldy'Z{y\y-fL)0{y'-fi) 



(53) 

Renormalizing in MS (generalized in the obvious way to cancel the 1/e 2 divergences), we 
find 



Z{y\y;P) 



2tt 



1 in* l N 

- + -l0g^-- 



mi 



a s (fl)C F 1 



71 



e (y' - y) 



S(y' -y) 
-0(y'-y)\ ■ 



(54) 



Note that the counterterm consists of a diagonal piece which is proportional to 5(y' — y), 
and an off-diagonal piece proportional to 8(y' — y). This latter terms mixes the operator 
0(y) with all operators 0(y') with y' > y. 
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Inserting the one-loop matrix element of the renormalized operator into (|45|) we find the 
the differential decay rate in the parton model b — > Xg'j 



dT 

dx 



T / dyC(y,x;fl){0(y;fl)) 



T J dyC(y,x;p,) 
a s C F 



—, — log — - 2 log — + . . . 

47T y m^ mi 



5(l-y + k, 



Air 



6(l-y + k, 



4 - 4 log 


2 




mf 




/lo 


+8 





1 -y + k + ) + 
log(l - y + k+)\ 



1-y + k, 



(55) 



One might worry about the appearance in (^) of logarithmic terms that depend on nib, since 
this scale has been integrated out and thus should not be present in the effective theory. 
These terms are due to our choice of factoring the heavy quark mass out of the soft scale 
m^l — y + k + ) by writing our expressions in terms of hatted quantities. The logarithms of 
mj cancel in the matching coefficient. 

The Wilson coefficients C(y,x;n) are determined by matching the collinear-soft theory 
onto LEET. In the collinear-soft theory, the Feynman diagrams for the forward scattering 
matrix element are shown in Fig. O. As with LEET, all divergences are regulated in 







+h.c. 



<c) (d) 
FIG. 11. Collinear-soft theory Feynman diagrams which contribute to the forward scattering 

amplitude through 0(a s ). 

dimensional regularization. Expanding the expression for the forward scattering amplitude 
obtained from these graphs in powers of (1 — x + k + ), we find for the differential decay rate 



dT 

dx 



Cl(fi)T < 



a s C F ( 2 P 2 , , i ft 
1 + — — log — + 5 log — + 
47r V m£ mi 



47T 



6(1 -x + A; H 



log(l — x + k^ 
1 — x + k+ 



5(1 - x + k^ 
1 



7 



x + k A 



(56) 
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Comparing Eqs. (|56| ) and (|55"D gives the short- distance coefficient C(y,x;fi). At tree level, 
the matching is trivial, and we write 



S(y-x) + ^C ( - 1 \y,x;P) 



C(y,x;p) = Cy{{j) 
where \x is the matching scale. At one loop, we find 



+ 0(a 2 s ), 



C^(y,x;P) = ( 21og 2 tL + 31og ^ + . . . ] S(y - x) - [ 3 + 41og^ 



r,2 



mi mi 



+A6(y - x) 
(2 log 2 
-4 



m b(y ~ x 

9{y-x) 

log 



/ iog(y-z) \ 

1,2 +3 log 



£ 2 



m 2 b {y-x) 



+ ... <%-x) 



y — a; m 2 (y — x) y — x 
At the scale p = m^y/y — x the logarithmic terms vanish, and we find 

C« (y, x; m fev ^^) = -3 + . . . . 

y - x 

The matching scale is therefore different for each operator 0(y). 



(57) 



P \ 0(y - x) 



(58) 



(59) 



dT 

dx 



B. Renormalization group 

The differential decay rate in LEET given in ( p5| ) may be written as 

a s C F 



T J dyC(y,x;p) < 
xS(l - y + k+) + 



1 - 



log^ 



P 2 



4vr \ mf (1 - y + &+) 2 
a s C F ( 4 



-2 log 



A 2 



47T 



log 



P 2 



ml{l - y + k + ) 2 , 
4 M 



1 — y + k + m 2 (l — y + k + ) 2 1 — y + k^ 



(60) 



and so the large logarithms in the matrix element of 0(y; p) vanish at the scale p = — 
y + k + ). (This expression looks highly singular, but as can be seen from (|5"5]), the delta 
functions combine with the other terms to form plus functions.) Thus, in order to sum all 
logarithms of /i we must continue to run the operator 0(y) in LEET. From (^) and 
we obtain the renormalization group equation 



d f 
^~rC(y, x;p) = / dy' 7(1/, y'\ p)C(y', x; p) , 
dfi J 

where j(y,y'; p) is the continuous anomalous dimension matrix 



(61) 



i{y,y';fr) 



a s (p)C f 
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P 2 



log— o - 1 <%' -y) 



mi 



(y' - y)- 



-0(y' - y) 



(62) 
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Solving ( pip analytically, however, is nontrivial and beyond the scope of this work ||. 
Instead, we may diagonalize the anomalous dimension matrix by taking high moments of 
the spectrum. This will allow us to compare our results to those of Refs. PJTOl. Note that in 
Refs. [§,|l(J] both leading and next-to-leading logarithms were resummed. This requires the 
two loop contribution to the 1/e 2 counterterm, the full one loop matching condition, and the 
two loop running of a s , none of which have been included here. As a result our calculation 
only resums the leading logarithms and a class of the subleading logarithms. However, it 
is straightforward to extract from the literature a resummation of exactly the same set of 
logarithms. 

To calculate the moments we set the residual momentum k to zero. (This residual 
momentum can easily be incorporated by boosting from the rest frame of the b quark, 
Pb = rribV, to the frame p = m b v + k). Taking moments unconvolutes the expression for the 
differential decay rate in LEET (f45|) and we obtain 



T(N) 











dx 


x N 


Jo 






= r 


f 


dx 


Jo 




= r 


f 


dz 


lo 





dx 



T dxx N_1 dyC(y - x; /i)(0(y; /i)) 

JO J-oc 



.JV-1, 



= T C(N;jl) (0{N;fl)) 



(63) 



where we have used 



C(y-x) = -C (l 

y \ 



X 



y , 



Q(y - x) 



(64) 



since C(y — x) just contains delta functions and plus distributions. Thus, the operator 
0(N; fi) is just a linear combination of the set of operators 0(y; fi). The matching from the 
collinear-soft theory onto LEET at tree level is trivial, and we find 



C{N-t) = C^) 



0(a 



(65) 



Determining C^(N; Hq) requires the one-loop expression of T(N) in the collinear-soft the- 
ory and the one- loop matrix element of 0(N; fi) between partonic states. The one- loop 
expression for the differential decay rate in the collinear-soft theory is given in d56|) . Setting 
k + to zero and taking moments we obtain 



r(iv) 




i 



21og ; 



N 
n 



VI N 1 2 A 2 p-i V 2 

7 log log — o - 5 log 



n Q 



+ 



(66) 



where no = e lE . This needs to be compared to the one- loop matrix element of (0(N; fj,)), 
which can be obtained from (|55|): 



(0(tf;/2))= / dyy N - l {0(y-fi)) 
Jo 



1 - 



47T 



,1 2 fi N ,i P-N 

4 log 2 — 4 log ■ 



m b n 



m b n 



+ ... . 



(67) 
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The one loop matching coefficient is now easily determined using (j63|) , (|66|) and (|67D and we 
find 



47T 



21og 2 ^ + 31og^— 



+ . . . . 



(68) 



At the matching scale fx = m b ^Jno/N all logarithms in this matching coefficient vanish. 
Furthermore, from fl6?D it is clear that the matrix element (0(N;fx)) contains no large 
logarithms of N at the scale fx = m b no/N. Thus we run in the collinear-soft theory from m b 



to rribyrio/N, perform the OPE, and run C(N; fx) from m b ^no/N to m b rio/N. 

The running of the coefficient Cy in the collinear-soft theory from the scale mj to the 



scale m b ^no/N is obtained by setting jl = m b ^n /N in The running in LEET is 

determined by the RGE for C (N; fx) 



lx-^C(N;fi)= 1 (N;fi)C(N;fi) 



where the anomalous dimension is given by 



l(N;fx) 



dz z N 1/ y(z; fx) 



a s (n)C F 



71 



1 - 2 log 



(— ) 



The solution to this equation is 



m b n . 



C{N ] ^) = C 2 v [m bx r-^) 



N 



n \ ( a s {m b ^) 



1£f( 1+ *z 21og^ 

00 \ ft"' 5 "0 



NJ \ a s (m b J§) 



N , 



(69) 



(70) 



(71) 



This sums perturbative logarithms of N into the coefficient C(N). We can then substitute 
this into (|63|) to obtain an expression for the resummed moments of the differential decay 
rate. 

Using the result for Cv{fx) given in (|39| ) and taking the matrix element of 0(N;fx) 
between hadronic states, we find find the resummed expression for large photon energy 
moments of the decay B —>■ X s 7 



T(N) = T f(N;m b n /N) 



'a s (m b -^\ 00 



( 5 ftX) 



O's 



a s (m b j^' 
a s (m by ^ 



2C F 



( l+TT 2 21og^- 



(72) 



Logarithms are explicitly summed in this expression and only long distance physics is con- 
tained in the function f(N;m b n /N). 

We can easily compare our results to those in the literature. A resummed expression for 
T(N) is given in Ref. @: 



r(N) = T f(N;m b /N)exp 



t — (2 —T c (ll) + T(m b y) + j(m b ^/y) 

Jn /N y \ Jm b y fX 



, (73) 
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where, at one loop, 

r e M = aMCF , m = - as{fi)Cp , 7 M = - ^.MCf (74) 

7T 7T ATX 

Note that the cusp anomalous dimension T c (/i) is the contribution to the anomalous dimen- 
sion from the 1/e 2 counterterm. Using only the one loop cusp anomalous dimension, tree 
level matching, and the one loop running of a s , Eq. (|73|) resums leading logarithms and 
the same class of next-to-leading logarithms we resum in our calculation. Performing the 
integrals in the exponent we reproduce ([72]) . Thus the approach presented here, based on an 
effective field theory, is in agreement with the factorization formalism approach for summing 
perturbative logarithms. 



V. CONCLUSIONS 

In the specific case of B — > X s 7 we have shown how Sudakov logarithms can be summed 
within an effective field theory framework. First we construct an intermediate theory which 
includes both collinear and soft degrees of freedom. By performing a one-loop calculation 
we show that this collinear-soft theory can be matched onto QCD at the scale m;, without 
introducing logarithmic terms into the short- distance coefficient. In addition we determine 
the one-loop anomalous dimension and solve the RGE. Next we integrate out collinear modes 
at the scale rriby/y — x by switching to LEET. We perform an OPE in powers of (y — x) 
which leads to the appearance of a nonlocal operator where two vertices are separated along 
the light-cone. The matrix element of this operator between B meson states is the structure 
function. We perform the OPE at one-loop in the collinear-soft theory and match onto the 
nonlocal operator in LEET. At the scale mby/y — x no logarithmic terms are introduced into 
the short- distance coefficient. 

In order to compare to the factorization formalism results in the literature we repeat 
our analysis for large moments of the decay rate. In this case we find that the collinear-soft 



theory matches onto LEET at the scale mbyn /N, and that there are no large logarithms in 
the matrix element of the bilocal operator at the scale mbUo/N. Using the renormalization 
group equations in the collinear-soft theory we sum logarithms of N between the scales mj 



and mbyrio/N. We then switch to LEET and sum logarithms of N between the scales 

m b\fno/N and mbUo/N. This sums all perturbative logarithms of N. We find that our 
result agrees with results presented in the literature. This gives us confidence that we have 
constructed the correct effective field theory. 

Though we have presented this work entirely in the context of B — > X s 7 our approach 
is general. It should be straightforward to apply the collinear-soft theory and LEET to 
other processes in which Sudakov logarithms arise. Furthermore, this approach could also 
be applied to exclusive decays, in which case one does not perform the final OPE onto 
LEET, but remains in the collinear-soft theory. This could be applied to recent results 
on factorization in nonleptonic decays ||32|| , as well as LEET-based relations between form- 



factors in decays to highly energetic final states [fffil . Since these latter results depend only 
on the spin symmetry of LEET, which is also present in the collinear-soft theory, they should 
remain valid in the present approach. 
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